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Q> ■ Abstract 



We discuss interacting fermion models in two dimensions, and, in particular, such 
that can be solved exactly by bosonization. One solvable model of this kind was 
proposed by Mattis as an effective description of fermions on a square lattice. We 
review recent work on a specific relation between a variant of Mattis' model and such 
a lattice fermion system, as well as the exact solution of this model. The background 
for this work includes well-established results for one-dimensional systems and the 
high- Tc problem. We also mention exactly solvable extensions of Mattis' model. 



^! 1 Introduction 

Q ■ The theory of metals developed in the 1950s-60s is the basis of an excellent qualitative 

■ description of many materials that exist in nature. This theory is based on the Fermi liquid 

picture of nearly-free quasiparticles and standard methods like diagrammatic perturbation 
theory; it is thoroughly discussed in most introductory textbooks on condensed matter 
physics (see e.g. Ref. pLj). 

ce- 
stui, this otherwise so successful theory is known to fail for systems with strong cor- 
relation effects, as is typically observed in low-dimensional systems. For example, in the 
case of fermions confined to one dimension a different picture has been developed that 
radically departs from that of a Fermi liquid. This picture is largely based on prototype 
models that can be solved exactly by analytical techniques, and which nevertheless allow 
to take into account certain types of correlation effects. One important such technique is 
known as bosonization, which amounts to rewriting interacting fermions in terms of free 
bosons; we will return to this in later sections. For fermions confined to two dimensions, 
an equally satisfactory description has not yet been achieved. In this review, we describe a 
particular approach that uses a combination of mean field theory and bosonization to study 
such fermion systems [21 [3l HI O [6l [7] . In particular, we discuss a class of two-dimensional 
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fermion models that can be solved exactly by bosonization, as well as their relation to mod- 
els describing fermions on a square lattice. We also discuss the so-called high-T^ problem 
providing one physical motivation for this work. 

Our plan is as follows. Section 11.11 introduces the models playing a central role in this 
paper, and Section [L2] contains a biased overview of the literature on bosonization in higher 
dimensions^ Section [2] contains a discussion of lattice fermion models and their relation 
to the high- Tc problem (Section 12.11) . As a pedagogical way of introducing the results and 
notation used in later parts, we also discuss the bosonization approach to one dimensional 
lattice fermion systems (Sections I2.2H2.3I) . Section [3] reviews results for a simple prototype 
two dimensional lattice fermion system obtained in Refs. [21 [3l HI [5] . Section H] describes 
extensions of the latter results to more complicated models [6|, [7]. The final Section [5] 
contains a somewhat subjective discussion of our results from a broader perspective. 



1.1 Exactly solvable QFT models 

The Luttinger model [HI El [lOl HI] is a prominent example of an exactly solvable quantum 
field theory (QFtI^) model that has become a prototype for correlated fermion systems in 
one spatial dimension (ID) [12j. This model describes two flavors of fermions, labeled by 
an index r = ± and by ID momenta k, and is defined by the Hamiltonian 

hl = J2 f ^^^^^^ + E / ^yrAp)M-p)JAp) (1) 

r=± r,r'=± 

with ipr^\k) fermion field operators defined by the usual canonical anti-commutator rela- 
tions (CAR), 

Jr{p) = J dk ■.^l{k)i,{k + p): (2) 

fermion density operators, Vryip) (suitabl^ two-body potentials, vp > the Fermi ve- 
locity, and the colons indicating normal ordering. The technique to solve this model is 
based on mathematical results collectively known as bosonization; see below. It 

is worth emphasizing that "exactly solvable" has a very strong meaning for this model: not 
only the eigenstates, but also all correlations functions can be computed analytically; see 
Ref. [13] and references therein. 

As first pointed out by Mattis [II], there exists a similar model describing fermions in 
2D and which is also exactly solvable by bosonization. This model describes four flavors 
of fermions that are labeled by two indices r, s = ± and by 2D momenta k = {ki, k2)', it is 

^The literature on higher dimensional bosonization is huge, and our discussion is far from exhaustive. 
^Note that by "QFT" we mean a quantum model with infinitely many degrees of freedom. 
3See e.g. [13]. 
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defined by the Hamiltonian 



r,s=i: 



r,s,r',s'=± 



/ 77r~T^'^r,s,r',s'{p)Jr.si—p)Jr',s'ip) 

J (27r)^ 



(3) 



with kg = (ki + 5^2) /"\/2, and Jr,s(p) = / (Pk -.i/jI g(k)ilj^ ^(k + p) : etc., similarly as above. 
Mattis, having the famous high- Tc problem of the cuprate superconductors [I5] in mind, 
proposed this model as an effective description of fermions on a square lattice. We believe 
that the Mattis model deserves to be better known: First, exactly solvable models in 2D 
are rare, and any physically motivated such example deserves to be studied from different 
points of view. Second, as proposed and elaborated by us in recent work [21 [Sj IH [5] , the 
exact solution of this model is a key part of a method to compute physical properties of 2D 
lattice fermion systems. Third, there exist exactly solvable extensions of the Mattis model 
that allow, for example, to study the effect of dynamical electromagnetic fields and spin on 
2D interacting fermion systems by exact solutions [6l [7] . 

A main part of this review is on the relation of the Mattis model and a model of fermions 
on a square lattice. This relation is based on a method to derive effective QFT models for 
such lattice fermion systems. This method, which we call (partial) continuum limit, is well- 
established in ID, and in the simplest case of spinless lattice fermions leads to the Luttinger 
model. The extension of this method to spinless fermions on a square lattice leads to a 
model that, in addition to the fermion fields ipr,}(k), r, s = ±, that describe the so-called 
nodal fermion degrees of freedom, also takes into account the so-called antinodal fermion 
degrees of freedom represented by operators i/jl^l(k), r = ±. While the energy dispersion 
relations of the nodal fermions are linear in the momenta, er.,s(k) = vprkg for r, s = ±, 
the antinodal fermions have hyperbolic dispersion relations, er.,o(k) = rcpk+k^ — /zq for 
r = ±, with computable constants cp and /io. The full Hamiltonian providing an effective 
description of 2D lattice fermions is the following extension of the Mattis model 

r,s=± ^ 

+ ^ d^k {cprk+k^ - /io) :?/'J^o(k)^r,o(k) : (4) 

r=± 

Vr,s,r',s' (P) Jr,s ( "P) Jr',s' (p) 

with two-body interaction potentials K-,s,r',s' determined by the underlying lattice fermion 
interactions [3]. The model defined by the Hamiltonian in (jlj) is a 2D analogue of the 
Luttinger model in the sense that it arises as a partial continuum limit of 2D lattice fermions, 
and that it is amenable to bosonization. We emphasize that this 2D Luttinger model is 
not exactly solvable; only the nodal fermion degrees of freedom can be mapped exactly to 
non-interacting bosons and thus be treated without approximations. Treating the antinodal 
fermions by mean field theory, one finds a significant parameter regime away from half filling 
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for which the antinodal fermions are gapped [1]. We proposed that, in this partially gapped 
phase, the low-energy physics of the system can be described by the Mattis model [4J. This 
motivated us to study the Mattis model from a mathematical point of view, and we found 
that the Mattis model is indeed exactly solvable in the same strong sense as the Luttinger 
model [S]. 

1.2 Bosonization in higher dimensions (review) 

The literature on bosonization in dimensions higher than one is quite extensive. Below 
we mainly discuss work whose bosonization methods lie closest to ours, while references to 
other approaches are briefly mentioned at the end. 

The first to apply bosonization to interacting, non-relativistic fermions in higher dimen- 
sion was Luther \l6i (see also Ref. |13j for contemporary work). For a discretized spherical 
Fermi surface in three dimensions, one can bosonize in the radial direction at each point 
on the surface, while treating the discretized angles as flavor indices [16j ; this is sometimes 
called radial- or tomographic bosonization. Using a generalized Kronig identity, one can 
then express the kinetic part of the Hamiltonian in terms of densities p!6l [T3] . Similarly, the 
fermion fields can be written in terms of charge shift operators (also called Klein factors), 
which depend on the radial direction, and an exponential of boson operators [13]. However, 
in this early approach, only density operators with momenta in the radial direction behave 
as bosons, and thus interactions with transverse momentum exchange cannot be treated 
exactly. 

Real interest in higher- dimensional bosonization came with the high- Tc problem and the 
prospect of finding non- Fermi-liquid behaviour in strongly correlated 2D models. In partic- 
ular, Anderson suggested that the 2D Hubbard model on a square lattice has a Luttinger 
liquid phase away from half-filling, and that this could be explored using bosonization meth- 
ods [17J. Whether or not Luttinger-liquid behaviour is possible in this model is still an open 
problem. Rigorous work on the renormalization group have shown that weakly-coupled 2D 
fermions with Fermi surfaces of non-zero curvature are in general Fermi liquids |18] (see 
also Ref. [19] and references therein). Most of the attention has therefore been on Fermi 
surfaces with "fiat parts", i.e. they contain portions that are straight (see for example 
Refs. [201 EI]). 

Prior to Anderson's suggestion, Mattis proposed a 2D model of fermions with density- 
density interactions that could be solved exactly using bosonization [T5j. The kinetic part 
contained four types of fermions with linearised tight-binding band relations on each side of 
a square Fermi surface (see ([3])), which is reasonable for a Hubbard- type model with short- 
range hopping and near half-filling. Unlike Ref. [16], the density operators were taken to 
be bosonic for all momentum exchange, although no details were provided. Rewriting the 
kinetic part in terms of densities, the Hamiltonian of the model could be diagonalized using 
a Bogoliubov transformation. 

In more recent work, Luther [22] also studied fermions on the square Fermi surface 
with linear tight-binding band relations. Let k\\ and k± denote the momenta parallel and 
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perpendicular to a face of the square Fermi surface. Following Ref. [16], one would treat k\\ 
as a flavor index, extend k± to plus and minus infinity, and introduce a Dirac sea in which 
all states k± < are filled. However, unlike the ID case, one should note that there is a huge 
degeneracy in choosing the accompanying flavor index to the unbounded momenta. One 
can for example do a Fourier transformation (change of basis) in the /c||-direction and then 
bosonize the fermions with a new flavor index x\\. These can be interpreted as coordinates 
for a collection of parallel chains aligned in the fc^-direction. In this way, Luther was able 
to include boson-like interactions with momentum exchange also in the A;||-direction [22] . 
The properties of this model were further investigated in Refs. [23l [2^ . 

Finally, we mention other work on rewriting interacting fermions as bosons, although 
these methods are somewhat different from that discussed here. Haldane ^25j formulated a 
phenomenological approach in which density fluctuations corresponding to momentum ex- 
change of states near an arbitrary D-dimensional Fermi surface are assumed bosonic. These 
ideas were further pursued in the works of Refs. [26l [271 [28] . Bosonization methods using 
functional integration and Hubbard- Stratonovich transformations have been developed in 
Refs. [29l[30]. There is also an approach based on non-commutative geometry [3T] . 



2 Motivation and background 

2.1 High- Tc problem and 2D lattice fermions 

The possible violation of Landau's Fermi liquid theory in models of strongly interacting 
fermions has been an actively researched problem for many years. Interest in this topic 
quickly grew with the discovery of high temperature superconductivity in the cuprates [TB] 
and the realization that these materials display many properties not described by Fermi 
liquid theory; see e.g. Ref. [32] for review. Early on, it was suggested that 2D lattice fermion 
models of Hubbard-type capture the strongly correlated physics of cuprates [33|, [311 ESI [36] . 
These models are appealing due to their apparent simplicity. For example, the much-studied 
2D Hubbard model can be defined by the Hamiltonian 

HtU = -tJ2Yl (4aS-,a + 4 + ^ ^^.t^i,!' ^i,o. = 4,aCi,a (5) 

a=t,i i 

with fermion operators 4 a labeled by sites i of a 2D lattice and a spin index a =t, i, 
and ^ sum over all nearest-neighbor (nn) pairs on this lattice (t > and U > are 

the usual Hubbard parameters). The Hubbard model is conceptually simple since, if one 
restricts to a lattice with a finite number of sites, the Hamiltonian in can be represented 
by a finite, albeit usually large, matrix^. 

However, it has proven to be very difficult to do reliable computations for 2D Hubbard- 
like models for lattice sizes and intermediate coupling values of interest for the cuprates 
(i.e. for U/t in the range between 2 and 10, say). Thus, despite of much work over many 

''For a lattice with Af sites, the size of this matrix is 4-^ x 4^. 
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years, no consensus has been reached on the physical properties of such models in large 
parts of the interesting parameter regime. To be more specific: one important parameter 
for Hubbard-like systems is filling u, which is defined as the average fermion number per 
lattice site. For the Hubbard model this parameter is in the range < u < 2. There is 
consensus that, at half-filling (i.e. for z/ = 1), the 2D Hubbard model describes a Mott 
insulator. 

The challenge is to study the 2D Hubbard model close to, but away from, half filling, 
which is a regime that is much less understood. Our approach has been to look for other 
models that describe the same low-energy physics but are more amenable to reliable com- 
putations. As we explain in Section 13. ![ our search for such models was motivated by mean 
field results for the 2D Hubbard model ^37j . Another important motivation and guide were 
ID lattice fermion systems that, different from 2D, are well understood also by numerical 
and analytical methods: In Sections I2.2H2.3I we review a particularly useful method in ID, 
namely to perform a suitable continuum limit to obtain a QFT model that can be stud- 
ied by bosonization. As already mentioned, this method is well-established in ID, but we 
present it in a way that makes our generalization to 2D [3J natural. For simplicity, much of 
our discussion in this paper is on a spinless variant of the Hubbard model defined by the 
Hamiltonian ^ 

Htv = -tYl + + Y ^ ^ ^^-^ 

with fermion operators c^^^ labeled only by lattice sites i; t > and > are model 
parameters. For this so-called t-V^model, the filling parameter is in the range < u < 1. 
Similarly as the Hubbard model, the 2D t-V^model is believed to describe an insulator at 
half-filling u = 1/2 and intermediate coupling values V/t > 0, and it has been a challenge to 
understand the model away from half-filling. We will also shortly describe generalizations 
of our results for the 2D t-l^model to the Hubbard model in Section 14.21 

2.2 ID lattice fermions and the Luttinger model 

We now describe a method to derive the Luttinger model from the ID variant of the so- 
called t-V^model defined in (E]). 

In ID, the lattice sites are i = 1,2, . . . ,M. Introducing a lattice constant a > and the 
"length of space" L = Ma, we define il){x) = a^^^'^Cj with the ID spatial positions x = ja. 
Due to the finite L, possible fermion momenta k are in the set A* = (27r/L)(Z +1/2) 
(we use anti-periodic boundary conditions). Using lattice Fourier transform, we express 
the Hamiltonian in ([6]) in terms of fermion operators ijj^'^^k) on the Brillouin zone of the 
lattic(i 

BZ=|A;gA*; -- < k < -] (7) 
I a a) 

^We use anti-periodic boundary conditions. 
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as follows 

Htv = [ dke{k)^\k)i){k) 

f\ f . .... (8) 

+ / dki--- dkiv{ki,- ■■ ,ki)i)\ki)i){k2)i)\k^)i){ki) 
Jbz Jbz 

with the band relation e{k) = —2tcos{ka) and the interaction vertex 

v{ki, . . . ,ki) = u{ki - /C2) ^ 5{ki - k2 + ks - k^ + {2n/ a)n) (9) 

where u{p) = aV cos{ap)/{2n)] we write Jgdk = Ylikes^'^'^ / ^) scaled Riemann sums 
(S* some subset of A*) and 5{ki — ^2) = [L/ {2'K)]5ki^k2 for scaled Kronecker deltas. Our 
normalization of the CAR is such that {'ip{ki),'ijj\k2)} = S{ki — ^2)- We note that these 
scalings are such that all formulas remain meaningful in the formal infinite volume limit 
L — )■ 00, but it is convenient to keep L finite during computations. We use a particle physics 
jargon and refer to a and L as UV- and IR- cutoff, respectively. 

Our aim is to modify the model so as to make it better amenable to analytical compu- 
tations but not (much) change the low-energy properties of the model. The strategy is to 
perform a continuum limit a — t- O"*", i.e., the lattice constant a is a small parameter, and 
sub-leading terms in a are assumed to be of lesser importance. However, this limit has 
to be taken with care: Physics folklore suggests that the low-energy physics of a fermion 
system is dominated by momentum states close to the Fermi surface, i.e., one should only 
change the Hamiltonian so that only momentum states are affected that are far away from 
the Fermi surface. To be more specific: In the non-interacting case = 0, the groundstate 
of the Hamiltonian in (jS]) is 

\n)= n V^^(^)lo) (10) 

with the chemical potential /i determined by the fermion density and |0) the vacuum state 
such that ■ilj{k)\0) = for all k. We now choose Q such that e{rQ/a) = /i for r = ±, 
i.e., the filled momentum states are in the range —Q/a < k < Q/a, and the two points 
±Q/a correspond to the Fermi surface. Note that filling in this situation is u = Q/ir, i.e., 
half-filling corresponds to Q = 7t/2. The groundstate \Q) is then fully characterized by the 
following conditions, 

ij±{±k)\n) = ■^i{Tk)\n) = yk>o (ii) 

where 'ipr{k) = '4^{Kr + k), r = ±, are fermion operators labeled my momenta close to the 
Fermi surface points Kr = rQ/a. 

To prepare for this limit it is convenient to rewrite the t-V^model as a model of two 
flavors of fermions. For that we divide the Brillouin zone in two subsets A* = {fc G 
BZ - Kr] r{k + Kr) > 0} such that 

[ dkf{k) = y^l dkf{Kr + k) (12) 
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for arbitrary functions / defined on the Brillouin zone. This shows that the fermion oper- 
ators iprik) are to be labeled by momenta k G A*, r = ±. Inserting this in (|8]) we obtain 
a Hamiltonian describing two flavors of fermion operators ipr{k), A; G A* and r = ±, with 
dispersion relations er{k) = t{K^. + k) and interaction vertices 

Vri,- ,ri{ki, . . . ,ki} = u{Kri — Kr^ + ki — ^2) 

X ^ 5{Kr^ - Kr^ + Kr., - + ki - k2 + ks - ki + {2n/a)n), (^3) 

nGZ 

i.e., Hfv = ^r=± /a* '^^ ^r{k)'ipl{k)'ii)r{k) + ■ ■ ■ (the dots indicate obvious interaction terms). 

We now modify the model as follows: First, we Taylor expand the dispersion relations 
er{k) = fi + rvpk + 0{a'^k'^), vp = 2tasm{Q), and ignore the higher order terms, i.e., we 
replace 

er{k) -)■ /i + rvpk. (14) 
Second, we replace the interaction vertex by 

Wri,-,r4(^l, • • • , A;4) -> u{Kr^ - Kr2)5{ki - k2 + k^ - ki) 

6{Kr, - Kr, + Kr, - Kr, + {2n /a)n), (^5) 

i.e., we expand u{Kr^ — K^^ + ki + k2) = u{{ri — r2)Q/a)[l + 0{a{ki — k2))] and ignore the 
higher order terms and, at the same time, modify allowed scattering terms in the interaction 
but only those where a{ki — k2 + k^ — k^) is "large". Since momenta k close to the Fermi 
surface are such that ka is "small", we expect that these changes are appropriate. 

Inserting these changes in the Hamiltonian one obtains, after some computations, 

Htv^y^f dkrvFktPl{k)Mk)+ [ ^ v{p)p^{p)p_{-p) (16) 



with 



Prip)= / dk^Plik)Mk + p) (17) 



au' 



vip) = gein - \p\a) (18) 

with g = 2aV sm'^{Q)/7i; here and in the following we use the notation A* = (27r/a)Z. We 
note that, in these computations, one can ignore additive constants and chemical potential 
terms, i.e. terms proportional to X]r=± /a* dk ipl{k)il)^{k) (since the former terms amount to 
constant shifts of the energy, and the latter to a renormalization of the chemical potential, 
which both are irrelevant). The final modification is to replace on the r.h.s. in f|T6|) and in 

m 

A* ^ A* (19) 

i.e., to remove the restriction on fermion momentum states by taking the limit a — )■ 0^ in 
the sets A*, r = ±. We note that it is this approximation which turns the model into a 



"9 below is the Heaviside function. 
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QFT: before we had a model with a finite number of fermion momentum states, but after 
(fT9|) we have a quantum model with an infinite number of degrees of freedom. One finds 
that, in order to obtain a model that is mathematically well-defined, one needs to normal 
order before making the change in (IT^ : defining 

■.A: = A - {Q,AQ) (20) 

for operators A and with the state \Q) fully characterized by f|TT]) . one can show that 

Mp)= [ dk :tljl{k)tijr{k+p): (21) 

J A* 

and 

HL = y^[ dkrvFk:^Pl{k)Mk):+ [ ^V{p)J+{p)L{-p) (22) 
^ J A* J A* 27r 

are mathematically well-defined. Note that normal ordering, before the replacement in 
f|T9|) . only amounts to dropping irrelevant additive constants and chemical potential terms. 

The Hamiltonian in fl22p makes mathematically precise a special case of the one in ([T]) 
(i.e. Vry = 5r-r'V /2). Note that the limit a — )■ 0''" was only partial and, in particular, it is 
important to keep a finite in the interaction potential in (1181) : taking a — >■ 0"*" amounts to 
V{p) — g (independent of p), which is exactly the limit where the Luttinger model becomes 
formally equal to the massless Thirring model |9]. However, as is well-known, this latter 
limit is delicate and requires a non-trivial multiplicative renormalization |38j . 

It is important to note that, in the derivation above, one assumes Q ^ -n j1. For Q = -n j1 
there are additional terms in the effective Hamiltonian that cannot be ignored and that, 
presumably, open a gap. This is consistent with Q = -r j1 corresponding to half filling. 



2.3 Bosonization in ID 

The exact solution of the Luttinger model relies on three mathematical facts: First, the 
normal ordered fermion densities do not commute (as one would naively expect) but obey 
the commutator relations 

\Jr{p)JrW\ = rKypl{p^pJ) (23) 

where the r.h.s. is an example of an anomaly. Second, the kinetic part of the Luttinger 
Hamiltonian can be expressed in term of the normal ordered fermion densities, 

/ dkrk:^\{k)^Xk):=T, I ^ :J,(p)X(-p):, (24) 

J A* J A* ^TT 

which is known as Kronig identity. Third, the fermion field operators can be computed as 

1 /"^/^ 



X exp ( ± / dp- e'P"-^IPl J±(p) ) e-'^^'rfx 
V A*\{0} P J 
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with unitary operators R± obeying R+R^ = —R^R+ and 



RrJr'{p)Rr = Jr'{p) — r5ry6pfi. (26) 

The commutator relations in (l23ll imply that the linear combinations 

fl(p) = i=(-J+(p) + J_(p)) 

"^l , , (P^O) (27) 

are boson operators with the usual canonical commutator relations (CCR) 

[n(p), ^\p')] = ~i6{p - p'), [U{p), &{p')] = Mp), ^\p')] = (28) 

with n^(p) = fl{—p) and similarly for $. Moreover, using fl2^ . one can express the Lut- 
tinger Hamiltonan solely in terms of normalized density operators, and inserting (127|) one 
finds 

Hl = ^ [ dp:(il- ^ip))n\pMp) + (1 + i{p))p'^\p)Hp)) ■■+■■■ (29) 
^ Ja*\{o} ^ ^ 

where the dots indicate zero mode terms that are 0(1/L) and that we suppress, to not 
clutter our presentation. The formula in fl29|) shows that the Luttinger Hamiltonian is 
equivalent to a free boson model that can be diagonalized by a boson Bogoliubov tran- 
formation. Moreover, since fl25l) expresses the fermion field operators in terms of boson 
operators, it is possible to compute all fermion correlation functions exactly by analytical 
methods; see e.g. Ref. [13] for details. 



3 Effective model for a 2D lattice fermion system 

In this section we review a series of papers where we propose [21 |3] and study [U [5] an 
effective model for a 2D lattice fermion system of Hubbard type. In these papers we treat a 
generalization of the 2D t-\^model where we also allow for a next-nearest neighbor hopping 
term proportional to t' but, for simplicity, we restrict ourselves here to the special case 
t' = 0. 



3.1 Mean field results in 2D 

Hartree-Fock theory is a variational method where the groundstate of the interacting system 
is approximated by the one of non-interacting fermions in an external potential, and this 
potential is determined such that the variational energy is minimized [39]. In mean field 
theory one further restricts to variational states that are translationally invariant, which 
often is adequate for translationally invariant models. For the t-l^ model, there are three 
mean field parameters, and one, which we call A, is proportional to the magnitude of the 
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density difference between neighboring sites, which is assumed to be constant A mean 
field solution with A > has the physical interpretation of a charge density wave (CDW) 
phase: the state is only invariant under translations by two sites, and it is usually insulating. 
A normal (N) phase corresponds to a state invariant under all lattice translations, i.e. A = 0, 
and such a state is typically metallic. It is important to note that mean field theory allows 
to compare the variational energy of three types of states: First, the pure N state, second, 
the pure CDW state, and third, a phase separated state where parts of the system are 
in the N and another part of the system in the CDW state |10l H]. If this third state is 
obtained it proves that the assumption underlying mean field theory is not correct: the 
true Hartree-Fock groundstate is not translationally invariant, and this is an indication of 
unusual physical properties that cannot be accounted for by mean field theory. 

Figure [Il^a) shows the mean field phase diagram for the 2D t-^model as a function of 
filling u and interaction strenght V/t, setting t = 1. As shown, the pure CDW state occurs, 
but only at strictly half filling, and the N state only in a region quite far away from half 
filling: in the shaded regions between mean field theory fails. We believe that this failure 
of mean field theory in a large region away from half filling is due to treating all fermion 
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(a) 2D t-V model 
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Figure 1: Comparison between the coupling {V) vs. filling (z/) mean field phase diagrams 
of the 2D t-^model (a) and the 2D Luttinger model (b) for t = 1. Shown are the charge- 
density-wave (CDW) and normal (N) phases as a function of coupling V and filling u at 
zero temperature. The regime marked by horizontal lines are mixed, i.e. neither the CDW 
nor the N phase is thermodynamically stable. The CDW phase in (a) exists only at half- 
filling u = 0.5. The CDW region in (b) corresponds to a partially gapped phase. The 
parameters used for (b) are k = 0.8 and Q/n = 0.45. (Figure reprinted from Ref. 4 with 
kind permission from Springer Science-f-Business Media.) 
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degrees of freedom in the same way. In the 2D Luttinger model these degrees of freedom are 
disentangled and thus can be treated by different methods. As shown in Fig. [TJ^b), mean 
field theory is applicable in a much larger part of the doping regime if used only for the 
antinodal fermions. As we proposed [U E], the CDW region in Fig. [U^b) corresponds to a 
partially gapped phase where the antinodal fermions are gapped and thus do not contribute 
to the low energy properties of the system, and this region can be described by the Mattis 
model. 



3.2 Partial continuum limit of the 2D t-Vm.ode\ 

We consider the model defined by the Hamiltonian in (E]), but now on a 2D diagonal 
square lattice with lattice constant a and size L, i.e., the lattice sites are x = (xi,X2) with 
Xi^2 integer multiples of a such that the possible fermion momenta are k = (A;i,/c2) with 
k± = (fci ± k2)^2 half-integer multiples of 21^ jL. We denote the set of all such fermion 
momenta as A*. 

The 2D Brillouin zone is a subset of a square (rather than of an interval in ID), 



BZ = <^kGA*; <A;i2<-L (30) 

and, similarly as in ID, the Hamiltonian in Fourier space is fully characterized by the 
dispersion relation 

e(k) = — 2t[cos(a/i;i) + cos(aA;2)] (31) 

and the interaction vertex 

t;(ki, . . . , k4) = u(ki - k2) ^ 5'(ki - k2 + kg - k4 + (27r/a)n) (32) 

neZ2 

with 'u(p) = a^\^[cos(api) + cos(ap2)]/(2vr)^, i.e., the Hamiltonian can be written as in ([8]) 
etc. but with the obvious replacements ?/'*^^)(fc) — )■ ^/''^'''^(k). 



etc. 



The non-interacting Fermi surface at half filling coincides with the tilted square given 
by the four line segments ki + sk2 = rn/{2a), r, s = ±, i.e., e(k) = on these lines. The 
midpoints of the sides of this square Fermi surface are {rQ,rsQ)/a with Q = 7r/2. We 
are interested in the model close to half-filling and assume that there is a Fermi surface 
containing four points 

Kr^s = {rQ,rsQ)/a, r, s = ± (33) 

close to these midpoints, i.e., Q ^ 7r/2 is a parameter. Close to these four points, the 
dispersion relation has a linear behavior, 

e{Kr,s + k) = -Atcos{Q) + rvFks + 0{\ak\'^), r, s = ± (34) 
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with vf = 2^/2sm{Q)ta. This suggests to approximate the dispersion relation in ( 131]) by a 
hnear one, keeping only the leading non-trivial term in This leads to a Fermi surface 
that includes line segments through the points K^^^, as indicated by dashed lines in Fig. [2l 
However, there are momentum states close to the non-interacting Fermi surface at half 
filling where this approximation is qualitatively wrong: the points K+ q = (^)0)/a and 
K_ o = (0,7r)/a on this Fermi surface are saddle points of the dispersion relation in f l3T|) . 
and 

e{Krfl + k) = -rcFikl-kl)/2 + 0{\ak\^), r = ± (35) 

with cp = 2ta^ . This hyperbolic behavior of the dispersion relation is important for the 
tendency of 2D lattice fermion system to develop a gap at half filling [H], and it certainly 
cannot be ignored. There are two more special points, namely = (tt, vr)/a and K_^2 = 
(0, 0) , corresponding to a minimum and maximum of the dispersion relation, respectively. 
However, the regions close to the latter points are far away from the assumed Fermi surface, 
and we expect that they therefore can be ignored. 
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Figure 2: Division of the Brillouin zone BZ in eight regions K^^s + A*^ marked as (r, s) for 
r = ±, s = 0, ±, 2. The eight dots mark the points K^^s, and the parameters are Q = 0.457r 
and K = 0.8. Note that ki^2 are only defined up to integer multiples of 27r/a, and thus the 
union of these regions cover the BZ exactly once. (Figure reprinted from Ref. 3 with kind 
permission from Springer Science-|-Business Media.) 



This suggests to us to write the 2D t-V^model as a model of eight flavors of fermion 
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operators V'r.slk) = + k), r = ± and s = 0, ±, 2, similarly as in ID. For that we use 

a partition of the Brillouin zone in eight regions K,. + A* ^ such that 

f d'kfik) = J2 E / ^"'^/(K.,s + k), (36) 

as indicated in Figure [2l As already mentioned, we refer to the fermions with s = ± as 
nodal and the ones with s = as antinodal. It is important to note that this partition 
depends on a parameter k in the range < k < 1: the antinodal regions are 

A*Q = A* = I k G A*; <k±< \ (37) 

I (X (XI 

with < K < 1, and this also determines the widths of the nodal regions: Ar,s=± only 
contains momenta k such that — vr/a < k^^ < T^/d with 

a = . (38) 

1 — K 

As discussed in more detail below, d is important since it serves as UV cutoff for the nodal 
fermions. 

We now modify the 2D t-V^model, written as a eight flavor model, similarly as in ID: 
First, the nodal- and antinodal dispersion relations are modified by ignoring the 0(|ak|")- 
terms in (n = 2) and [n = 3), respectively. Second, the interaction vertices are 
modified by using the obvious generalization of (IT^ and (ITB]) to 2D. Third, the restriction 
on the nodal momenta orthogonal to the assumed Fermi surface is removed, i.e., 

A;,^±^A: = (kG A*; -^<A;_, (39) 
I. a a J 

with d in fl38|) . and, at the same time, all terms in the Hamiltonian involving fermion flavors 
with s = 2 are ignored. Note that this third change only modifies fermion degrees far away 
from the assumed Fermi surface. Similarly as in ID, it is the change in f l5I?]l that leads to 
a QFT model, and this model is mathematically well-defined if this change is done after 
normal ordering [3]. By lengthy computations one finds a Hamiltonian which, remarkably, 
only contains interaction terms involving densities 

:p,,,(p):=/ d^k.f d% :V'L(ki)V'..(k2) :5'(ki-k2 + p) (40) 

for r = ±, s = 0, ±. A final change is needed to obtain the 2D Luttinger model: for reasons 
explained below, we replace in the Hamiltonian 



Pr,s=± 



;p):^X,,(p)^ /" d'k^ [ d'k2 ■^iMi^^.^k,) : 

5^{ki - k2 + p + 2nne_s/a) 



X 

nez 
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with e± = (l,±l)/v2. Note that this last change amounts to adding umklapp terms to 
the nodal densities. 

The changes described above lead to [3] 

HtV ^£o + Hn + Ha + Hna (42) 

with the nodal part of the Hamiltonian 



r,s=±-^ "'^* r,sy,s'=± 
Xg {6s,s'Sr-r' + Ss-s') Xs (p)Xs' (p) -^r-,s ("P) (p) > 

the antinodal part 

Ha = ^ d^k {-rcpk+k^ - /io) : V'J,o(k)V'r,o(k) : 
— L-JA* 



72 



dp 

+ / 7T^9 ■■ P+,o(-p) P-,o(p) 

A* [^^) 



(43) 



(44) 



and the nodal-antinodal interactions 

Hna= J^-^ 22 9Xs{p) -Prfii-P)- Jt'Ap), (45) 

with the coupling constants g = 2Va? sin^(Q), g = 2Va? , and the cutoff functions Xs(p) = 
Qiixja — |p_s|)6'(K7r/(-\/2a) — \ps\)\ see Refs. [3l H] for further details, including explicit for- 
mulas for the constants and /io, and how filling v depends on Q and k. The Hamiltonian 
H = Hn + Ha + makes mathematically precise a special case of the one in (jlj). 

Similarly in ID, one assumes Q ^ tx/2 in the above derivation, and for Q = 7r/2 
additional nodal interaction terms appear that cannot be bosonized in a simple manner 
and that, presumably, open a nodal gap. This is consistent with Q = 7i/2 corresponding to 
half filling. 



3.3 Bosonization in 2D 



The 2D analogues of the identities in 



and dMD are [S] 

2vrp, 



d^krk, -.tfjl {k)ilj (k): = dn 

At J A [2TT 



(5^(p + p' - 2nne_Jd) 

(i^p 



■Jr,si-p)Jr,siP)- 



(46) 



(47) 



with r,s,r's,'= ±; see Ref. Proposition 2.1 for a mathematically precise formulation. 
Moreover, there exists a generalization of (12^ to 2D; see Ref. [3], Proposition 2.7. 
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a 1 



Similarly as in ID, the first of these identities implies that the operators 

n.(p) = [-J+Ap) + J-Ap 
•^^(p) = V J~~ (j+Ap) + J-Ap 

are boson operators with the usual CCR, and using these operators, the nodal- and mixed 
parts of the 2D Luttinger Hamiltonian can be written as [3] 

H- = Y^L ' ' ^'^^^^ nt(p)n,(p) + (1 + 7.(P)) P^$I(P)*.(P) 

^' ^ (49) 

+7(p)p+P-$I(p)$,(p) 1 : + ... 



and 

72^ 



/ (Iti^ ^ yY^Ps9Xsip) -prA-p)- ^sip), (50) 

with 7,(p) = 7Xs(p), 7(p) = 7X+(p)X-(p), and 7 = 1/(1 - k) sin(Q)/(27rt); as before, the 
dots indicate zero mode terms that are 0{1/L). 



3.4 Relation to the Mattis model 



As described in the previous section, the 2D Luttinger Hamiltonian can be mapped exactly 
to a Hamiltonian describing non-interacting bosons coupled linearly to interacting antinodal 
fermions. It is possible to integrate out these nodal bosons exactly and thus obtain an 
effective model for the antinodal fermions only. In Ref. [4] we studied this latter model 
by conventional approximation methods including mean field theory, and we found a non- 
trivial region away from half filling where the antinodal fermions are gapped. We assume 
that the gapped fermions do not contribute to the low energy properties of the system, and 
we therefore proposed that the Mattis model describes the low energy physics of the 2D 
Luttinger model in this partially gapped phase. 



3.5 Solution of the Mattis model 



Since the nodal part of the 2D Luttinger Hamiltonian is quadratic in boson operators, it 
can be diagonalized by a boson Bogoliubov transformation. From this, we can compute 
the ground state- and free energy of the nodal (Mattis) model [5]. Furthermore, there is 
a natural generalization of the boson-fermion correspondence formula in fl25|) to the two- 
dimensional case |5]. This formula allows to compute all fermion correlation functions 
exactly using standard results for free bosons. For example, we found that the fermion 
two-point function has algebraic decay for intermediate length scales, and with non-trivial 
exponents that depend on the coupling constants. 
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The result for the two-point function may be interpreted as a signature of Luttinger- 
hquid behavior in the Mattis model. We emphasize however that this does not automatically 
imply Luttinger-liquid behavior in the 2D t-V^model of lattice fermions, even if the antinodal 
fermions are gapped. For this to hold true, one first needs to investigate the approximations 
introduced when deriving the 2D Luttinger model from the lattice system. 

4 Extensions of the Mattis model 

There are several ways to obtain new exactly solvable models of interacting fermions by 
extending the Mattis model. Below we shortly describe two such extensions [EIC]- 

4.1 A 2+ ID quantum gauge theory with interacting fermions 

In Ref. |6], we consider the quantum gauge theory model obtained by minimally coupling 
the fermions in the Mattis model to a two-dimensional dynamical electromagnetic field. 
The gauged model shares many features with the so-called Schwinger model |42J of (1+1)D 
quantum electrodynamics with massless Dirac fermions; see also Refs. [^Bl] . The anomaly 
in ( H6|l and the requirement of gauge invariance on the quantum level leads to a bare mass 
term proportional to the electric charge for the electromagnetic field, similarly as for the 
Schwinger model. 

The Hamiltonian of the gauged model can be bosonized and subsequently diagonalized; 
one finds two gapped, or massive, boson modes and one gapless mode. We also computed 
the linear response of the magnetic field to an external current and found that there is a 
Meissner effect. 

4.2 Partial continuum Hmit of the 2D Hubbard model 

Our work on the 2D t-V^model of spinless fermions can be extended to the Hubbard model 
(see (^). In Ref. [7j, we applied the partial continuum limit to the 2D Hubbard model 
on a square lattice. This again leads to an effective QFT model of nodal fermions coupled 
to antinodal fermions. The nodal fermions can be bosonized using either abelian- or non- 
abelian |45j methods, and in the latter case we obtain a natural 2D analogue of a Wess- 
Zumino-Witten model. Furthermore, by a specific truncation of the nodal part of the 
effective Hamiltonian, we obtain a spinfull variant of the Mattis model that is exactly 
solvable. The fundamental excitations of this model separate into independent spin- and 
charge degrees of freedom. 
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5 Discussion 



The high- Tc problem has been an outstanding challenge in theoretical physics for many 
years; see for example Ref. |16] for recent commentaries by leading experts on the status of 
this field. One important hypothesis underlying much of the work on the high- Tc problem 
is that the 2D Hubbard model is an adequate prototype model for the cuprates. From a 
theorist's point of view this hypothesis is highly appealing due to the (deceptively) simple 
and aesthetic form of the model. Furthermore, the model is considered "adequate" in the 
sense that, (i) it is possible to derive it using physical arguments from a more realistic 
description and, (ii) it contains enough physics to capture the qualitative properties of the 
cuprates, while better quantitative agreement can be obtained by straightforward extensions 
of the model. Still, while there seem to be little doubts about (i) and (ii), these criteria are 
not enough; an adequate model also requires that, (Hi) it should be amenable to reliable 
computations and thus allow for experimental predictions. This criterion, which is obviously 
the most important one, remains to be fulfilled. 

The results reviewed in this paper are part of a program aiming at finding a better 
balance between (i)-(iii) by modifying a 2D Hubbard-like lattice fermion model. Our 
emphasis has been on modifications that lead to exactly solvable models since, for such 
models, (Hi) is always true. However, it is important to note that, even though we have 
proposed a derivation of such a model from a Hubbard-like model, we cannot claim to have 
shown that (i) is fulfilled: our derivation contains steps that need to be substantiated by 
mathematical arguments (these steps are spelled out in Refs. [3] and [5]). We also stress 
that, up to now, (ii) remains open: we have not systematically investigated whether the 
Mattis model, or its various extensions, can describe the cuprates. Much remains to be 
done in our program. 

Acknowledgements 

This work was supported by the Goran Gustafsson Foundation and the Swedish Research 
Council (VR) under contract no. 621-2010-3708. We thank Farrokh Atai for carefully 
reading the manuscript. 

References 

[1] G. D. Mahan, Many particle physics (Plenum Press, 1981). 

[2] E. Langmann, "A two dimensional analogue of the Luttinger model" , Lett. Math. Phys. 
92, 109 (2010). 

[3] E. Langmann, "A 2D Luttinger model", J. Stat. Phys. 141, 17 (2010). 

[4] J. de Woul and E. Langmann, "Partially gapped fermions in 2D", J. Stat. Phys. 139, 
1033 (2010). 



18 



[5] J. de Woul and E. Langmann, "Exact solution of a 2D interacting fermion model", 
larXiv:1011.1401 (to appear in Comm. Math. Phys.). 

[6] J. de Woul and E. Langmann, "Gauge invariance, correlated fermions, and Meissner 
effect in 2+1 dimensions" . larXiv: 1107708911 

[7] J. de Woul and E. Langmann, "Partial continuum limit of the 2D Hubbard model", to 
appear. 

[8] S. Tomonaga, "Remarks on Block's method of sound waves applied to many-fermion 
problems". Prog. Theor. Phys. 5, 544 (1950). 

[9] W. Thirring, "A soluble relativistic field theory", Ann. Phys. 3, 91 (1958). 

[10] J. M. Luttinger, "An exactly soluble model of a many-fermion system", J. Math. Phys. 
4, 1154 (1963). 

[11] D. C. Mattis and E. H. Lieb, "Exact solution of a many-fermion system and its asso- 
ciated boson field", J. Math. Phys. 6, 304 (1965). 

[12] F. D. M. Haldane, ""Luttinger liquid theory" of one- dimensional quantum fluids: L 
Properties of the Luttinger model and their extension to the general ID interacting 
spinless Fermi gas", J. Phys. C 14, 2585 (1981). 

[13] R. Heidenreich, R. Seller, and D. A. Uhlenbrock, "The Luttinger model", J. Stat. Phys 
22, 27 (1980). 

[14] D. C. Mattis, "Implications of infrared instability in a two-dimensional electron gas", 
Phys. Rev. B 36, 745 (1987). 

[15] J. G. Bednorz and K. A. Miiller, "Possible highTc superconductivity in the BaLaCuO 
system", Z. Physzk B 64, 189 (1986). 

[16] A. Luther, "Tomonaga fermions and the Dirac equation in three dimensions", Phys. 
Rev. B 19, 320 (1979). 

[17] P. W. Anderson, "Luttinger-liquid behavior of the normal metallic state of the 2D 
Hubbard model", Phys. Rev. Lett. 64, 1839 (1990). 

[18] M. Disertori and V. Rivasseau, "Rigorous proof of Fermi liquid behavior for Jellium 
two-dimensional interacting fermions", Phys. Rev. Lett. 85, 361(2000). 

[19] M. Salmhofer, Renormalization: An introduction (Springer, New York, 1999). 

[20] A. T. Zheleznyak, V. M. Yakovenko and I. E. Dzyaloshinskii, "Parquet solution for a 
fiat Fermi surface", Phys. Rev. B 55, 3200 (1997). 

[21] V. Mastropietro, "Luttinger liquid fixed point for a two-dimensional flat Fermi surface" , 
Phys. Rev. B 77, 195106 (2008). 

[22] A. Luther, "Interacting electrons on a square Fermi surface", Phys. Rev. 5 50, 11446 
(1994). 



19 



[23] J. O. Fjaerestad, A. Sudbo and A. Luther, "Correlation functions for a two dimensional 
electron system with bosonic interactions and a square Fermi surface", Phys. Rev. B 
60, 13361 (1999). 

[24] O. F. Syljuasen and A. Luther, "Adjacent face scattering and stability of the square 
Fermi surface", Phys. Rev. B 72, 165105 (2005). 

[25] F. D. M. Haldane, "Luttinger's Theorem and bosonization of the Fermi Surface", 
in Perspectives in Many-Particle Physics — Proceedings of the International School of 
Physics Enrico Fermi, Course CXXI, ed. R. Broglia and J. R. Schrieffer (North Holland, 
1994) p. 530. 

[26] A. H. Castro Ncto and E. Fradkin, "Bosonization of the low energy excitations of Fermi 
liquids", Phys. Rev. Lett. 72, 1393 (1994). 

[27] A. Houghton and J. B. Marston, "Bosonization and fermion liquids in dimensions 
greater than one", Phys. Rev. B 48, 7790 (1993). 

[28] R. Hlubina, "Luttinger liquid in a solvable two-dimensional model", Phys. Rev. B 50, 
8252 (1994). 

[29] P. Kopietz and K. Schonhammer, "Bosonization of interacting fermions in arbitrary 
dimensions", Z. Phys. B 100, 259 (1996). 

[30] J. Frohlich, R. Gotschmann and P. A. Marchctti, "Bosonization of fermi systems in 
arbitrary dimension in terms of gauge forms", J. Phys. A 28, 1169 (1995). 

[31] A. P. Polychronakos, "Bosonization in higher dimensions via noncommutative field 
theory", Phys. Rev. Lett. 96, 186401 (2006). 

[32] D. A. Bonn, "Are high-temperature superconductors exotic?". Nature Physics 2, 159 
(2006). 

[33] P. W. Anderson, "The resonating valence bond state in La2Cu04 and superconductiv- 
ity". Science 235 , 1196 (1987). 

[34] V. J. Emery, "Theory of high- Tc superconductivity in oxides", Phys. Rev. Lett. 58, 
2794 (1987). 

[35] C. M. Varma, S. Schmitt-Rink, and E. Abrahams, "Charge transfer excitations and 
superconductivity in ionic metals". Solid State Comm. 62, 681 (1987). 

[36] F. C. Zhang and T. M. Rice, "Effective Hamiltonian for the superconducting Cu ox- 
ides", Phys. Rev. B 37, 3759 (1988). 

[37] E. Langmann and M. Wallin, "Mean field magnetic phase diagrams for the two dimen- 
sional t-t' -U Hubbard model", J. Stat. Phys. 127, 825 (2007). 

[38] K.G. Wilson, "Operator-product expansions and anomalous dimensions in the Thirring 
model", Phys. Rev. D 2, 1473 (1970). 



20 



[39] V. Bach, E. H. Lieb, and J. P. Solovej, "Generalized Hartree-Fock theory and the 
Hubbard model", J. Stat. Phys. 76, 3 (1994) 

[40] E. Langmann and M. Wallin, "Mean-field approach to antiferromagnetic domains in 
the doped Hubbard model", Phys. Rev. B 55, 9439 (1997). 

[41] H. J. Schulz, "Fermi-surface instabilities of a generahzed two-dimensional Hubbard 
model", Phys. Rev. B 39, 2940 (1989). 

[42] J. Schwinger, "Gauge invariance and mass. H", Phys. Rev. 128, 2425 (1962). 

[43] N. S. Manton, "The Schwinger model and its axial anomaly", Ann. Phys. (N. Y.) 159, 
220 (1985). 

[44] H. Grosse, E. Langmann and E. Raschhofer, "The Luttinger-Schwinger model". Annals 
of Phys. (N.Y.) 253, 310 (1997). 

[45] E. Witten, "Non-abelian bosonization in two dimensions". Comm. Math. Phys. 92, 
455 (1984). 

[46] J. Zaanen et al, "Towards a complete theory of high Tc" , Nature Physics 2, 138 (2006). 



21 



